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THE  GRAVITY-CAPILLARY  WAVE  INTERACTION 
APPLIED  TO  WIND-GENERATED,  SHORT-GRAVITY  WAVES 


1 


! 


INTRODUCTION 

The  second-order,  gravity-capillary  wave  interaction  elucidated  by  Valenzuela  and 
Laing  [1]  has  recently  been  found  to  play  a major  role  in  the  equilibration  of  short-gravity 
waves  in  a wind-wave  tank  [2] . The  exact  equations  for  the  interaction  require  the  numer- 
ical evaluation  of  several  integrals  and  thus  do  not  lend  themselves  easily  to  the  development 
of  a simple,  intuitive  model  of  the  interaction.  Furthermore,  the  exact  theory  is  developed 
for  irrotational  waves  and  includes  no  wind-speed  dependence  except  in  the  form  of  the 
directional-wave  spectrum. 

In  this  report,  the  gravity -capillary  wave  interaction  is  reexamined  for  the  special  case 
where  one  component  of  the  resonant  triad  is  in  the  short-gravity-wave  range.  Wind-speed 
effects  are  included  in  an  ad  hoc  manner  by  forcing  the  dispersion  relation  to  conform  to 
that  found  in  microwave  scattering  experiments  [3] . When  ak"4  capillary  wave  spectrum  is 
used,  the  result  is  a simple,  analytical  expression  for  the  energy  transfer  which  agrees  well 
with  numerical  calculations  in  the  no-wind  case.  The  expression  shows  that  the  interaction 
may  be  envisioned  as  one  between  a short-gravity  wave  and  a capillary  wave  whose  group 
speed  equals  the  phase  speed  of  the  gravity  wave.  The  interaction  dies  out  for  long-gravity 
waves,  having  a fc5/2  dependence  in  the  no-wind  case. 


GENERAL  THEORY 


Valenzuela  and  Laing  developed  the  theory  from  the  well-known  equations  for  velocity 
potential  and  surface  deviation  in  which  they  included  surface  tension.  These  two  functions 
were  expanded  to  third  order  in  wave  slope,  and  each  order  was  considered  a Fourier  series. 
From  the  solutions  of  differential  equations  for  different  orders,  they  computed  the  energy 
transfer  up  to  fourth-order  terms.  Their  results,  corrected  by  a factor  of  two  in  the  energy 
transfer,  are  (symbol  definitions  follow): 
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where  #x  and  #2  give  the  signs  associated  with  cox  and  w2  in  the  resonance  condition 

w2  ” = w3  » k2  ^ kl  = k3  • 

The  angle  0 between  kx  and  k3  is  given  by 
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where  upper  and  lower  signs  in  Eq.  (4)  go  with  those  in  Eq.  (3). 

The  following  additional  definitions  apply  in  the  above  equations: 
k(  is  the  wave  number  of  ith  wave, 
o>i  m y/gkj  + Tfe-3  is  the  angular  frequency  of  ith  wave, 

Oy  is  the  angle  of  ith  wave  with  symmetry  axis, 

-S(Aty,«y)  is  the  wave-number  spectrum  of  ith  wave  in  cm4.  Multiplied  by 
pw«y/kj  this  gives  the  energy  spectrum  (pw  = density  of  water), 

g ■ 980  cm/s2, 

T - 74  cm3  /s2, 

km  = s/gfT  - 3.64  cm"1. 

The  sum  over;  in  Eq.  (1)  accounts  for  the  ambiguity  in  Oj  and  a2  due  to  the  sign  ambiguity 
in  0 resulting  from  Eq.  (4). 


WIND-DRIVEN,  SHORT-GRAVITY  WAVE  THEORY 

We  now  apply  these  results  to  the  case  in  which  the  subscript  3 refers  to  a short-gravity 
wave.  In  this  case,  only  difference  interactions  exist;  i.e.,  Eq.  (3)  can  be  satisfied  only  for  the 
lower  signs.  Furthermore,  to  include  wind-speed  effects  in  the  dispersion  relation,  we  write 

Wy  = v/gfey  + Tfe3  + ky-Uy  = WyQ  + fcy*U,  (5) 

where,  following  Plant  and  Wright  [3] , 


. , . rp~  u*  / 0.277 


Here  pa  and  are  the  densities  of  air  and  water,  Ug  is  the  surface  water  velocity  at  an  air 
friction  velocity  u, , and  z0  is  the  roughness  length  at  u, . These  parameters  have  all  been 
determined  experimentally. 

The  assumption  that  k3  refers  to  a short-gravity  wave  requires  that  kx  and  k2  refer  to 
the  capillary-wave  range.  Thus, 

gk 3 > , gkx  < Tk\  , and  gk2  < Tk\  . 


(7) 
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This  implies  that  ko/k1  = k3/k,  is  small  and  that  <*)30/w10  = y/ic^/K^  . Thus  («30/<o10)2 
<c3/«1 . In  the  following  development,  we  shall  drop  terms  which  are  higher  than  first  order 

Ko/K,. 


3'  1 ' 


The  resonance  condition,  Eq.  (3),  may  be  written 


*2  ' kl  + k3> 
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The  geometry  of  this  interaction  is  shown  in  Fig.  1.  Equations  (7)  and  (8)  lead  immediately 
to  the  following  relations: 
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Fig.  1 — Geometry  of  the  gravity -capillary  wave 
interaction  applied  to  a abort -gravity  wave 
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where  we  have  used  the  fact  that  = l/2 . Solving  for  cos  0,  we  obtain 

„ 2*ilc30  + (U3 -l/jjcosaa) 

“S(J 53^ (ID 

where  c3q  = \/glk3  is  the  gravity-wave  phase  speed  in  the  absence  of  wind.  Equation  (11), 
then,  replaces  Eq.  (4)  when  k3  refers  to  a short-gravity  wave.  Note  that  the  sign  ambiguity 
still  exists  for  0. 

We  must  now  simplify  Eq.  (1).  First  note  that  difference  interactions  require  that 
*1  * ~ “d  *2  = +-  Thus  = 0 and  T ^ = 0,  so  that  only  the  first  integral  in  Eq.  (1) 

contributes.  Let  us  concentrate  first  on  D k"  • *k  which  we  shall  simply  call  D.  Substituting 
Eqs.  (3)  and  (9)  into  Eq.  (2)  and  using  Eq.  ^7)  yields 

D • 2 2wS*l(*l  + *3co®0)  “ w3*l(fel  + k3c°sp)  + 


<o|  I k\kzcos&  + (*!  /t3cosP)fc1fc3co8g|  1 
*3  [ W1  Wj  + w3  j|  I 


^ = + »<»>3fciC08fl  + + k3COS0 

2 l Wj  Wj  + u>3 


Thus,  to  first  order  in  Kj/Xj  , 
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The  expression  for  the  energy  transfer  may  be  written  to  this  same  order  as 
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(16) 


Since  all  terms  inside  the  braces  are  of  order  k3/k^  or  higher,  we  need  to  evaluate  ) 
only  to  lowest  order: 


_ 47TQ}1kfcOS2p 
1 ~ 3lsin0l 


3co10  + 3fe1t/1cosa'j 


,(/) 


3cj10  + 2fe1f/1cosa0) 


, Icos0l  < 1.  (17) 


When  Icos0l  > 1,  ) = 0.  Thus,  the  lower  limit  of  the  integral  in  Eq.  (16)  may  be  replaced 

by  fej,  the  value  of  kx  when  cosfi  = 1.  From  Eq.  (11), 

3 y/Tkt  /2  + l/j cosag  = y/g/k3  + U3 cosa3.  (18) 


The  wind-speed  dependence  of  the  term  in  brackets  in  Eq.  (17)  is  very  weak.  For  conditions 
of  practical  interest,  it  varies  from  about  1 to  1.3.  Thus  we  shall  evaluate  it  at  fe(  for 
aj  = 0 and  denote  it  as  A . Then. 
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EVALUATION  FOR  k * CAPILLARY-WAVE  SPECTRA 


Capillary-wave  spectra  measured  in  wind-wave  tanks  fall  off  approximately  as  fc“4  and 


have  an  angular  dependence  close  to  cos2a  [4] . Thus  we  now  assume  that 
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for  the  capillary  waves.  Then,  if  Eq.  (11)  is  used  to  evaluate  Wg/cjj , Eq.  (16)  yields 
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Fig.  1 shows  that  = a3  ± (3  which  accounts  for  the  sum  over  j in  the  above  equations. 


Two  cases  are  of  particular  importance,  namely,  when  a3  = 0,  and  when  an  integral  over 
a3  is  performed.  For  a3  = 0,  we  note  that  a1  = ± /3  so  that  the  sum  over  ; causes  the  term 
involving  sin  aW  to  drop  out.  Then  we  have 


dS(k  3,0)  (m  16AirC2coxk3cos&pdk1 


dt  fk  3fe6 1 sin/31 


(23) 


S(k3, 0) 


i: 


“ 56i4jrCcu3fe3cos5/3  dkx  f°°  8AirCoJ3k3U1cos4(}(2  - 3cos2P)dkj 


3k2  I sin/31 


L 


3u)1ki  I sin/3 1 


Note  that  Eqs.  (11)  and  (18)  indicate  that 

cos/3  = y/k//k1  , and  sin/3  = y/l~  kl/k1 


(24) 


If  we  use  this  and  let  Ul  = Uf  which  is  independent  of  fc1 , Eq.  (23)  becomes 
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Now,  the  only  troublesome  integral  is  the  last  one  wiiich  contains  in  the  denominator. 

Using  Eq.  (24),  we  may  write 

o>ilk-5'2  = k-*[T112  +kjl2Uskl1]-1  (26) 

As  fej  goes  from  k{  to  °°,  the  term  in  brackets  goes  from  [T112  + 1 to  T 1/2  while  k 4 

goes  rapidly  from  fe(~4  to  zero.  Thus,  the  main  contribution  of  the  term  in  brackets  will  be 
near  kt,  so  we  approximate  cjj  by 

cox  i k\12  [T112  + t/gfe,“1/2].  (27) 

Once  again  note  that  no  approximation  is  involved  if  the  wind  speed  is  zero. 

The  integrals  in  Eq.  (25  )are  now  easily  evaluated  to  give 
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If  the  wind  speed  is  zero,  this  equation  becomes 


dS(fc3,0) 


= 205C2rV7/2fe2/2  - i29crg"1/2fc6/2s(fe3,0) 


This  expression  has  been  compared  with  the  original  theory  by  using  a computer  program 
developed  by  Valenzuela  and  Laing  [1] ; we  used  S(k3, 0)  = 0.0064fe34  for  k3  > 0.364  cm-1 
and  zero  otherwise.  The  results  are  given  in  Table  1. 
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Table  1 — Comparison  of  Energy  Fluxes  Predicted  by  Original  Theory  and  by  Eq.  (29) 


k3 

(cm'1 ) 


9S(fe3,0)/9t  (Eq.  1) 

9S(fc3,0)/9t  (Eq.  29) 

(cm4) 

(cm4) 

-4.86  x i<r2 

-5.71  X 10'2 

-1.60  X 10'2 

-2.02  X 10"2 

-8.09  X 10'3 

-1.10  X 10'2 

-4  96  X 10'3 

-7.10  X 10'3 

-3.40  X 10'3 

-5.10  X 10'3 

Now  let  us  consider  the  second  case  mentioned  above,  that  of  integrating  the  energy 
transfer  given  by  Eq.  (22)  over  a3 . We  assume  the  following  form  for  S(k3  ,a3 ) in  conformity 
with  wave-tank  experiments: 


S(k3,a3) 


S(k3)cos2a3 

0 


lofgl  < jt/2  , 
lcr3 1 > w/2  . 


Furthermore,  in  order  to  simplify  the  integrations  over  a3 , we  assume  that  a(  = 0 in  all  dis- 
persion relations  so  that  all  dependence  on  a3  is  in  the  spectral  form.  In  view  of  the  ad  hoc 
manner  in  which  wind  speed  is  being  incorporated  in  the  theory,  this  approximation  is 
probably  satisfactory.  As  usual,  the  approximation  becomes  exact  for  zero  wind  speed. 

The  integrals  over  a3  are  now  straightforward  if  one  recalls  that  aj  = c*3  + 0 and 
a2  = a3  — 0*  We  obtain 
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Then,  letting  U.  = Uf  and  using  Eq.  (26)  for  in  the  last  integral  and  Eq.  (24)  for 
an  p and  cos  p,  we  obtain 
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Finally,  evaluating  the  integrals  and  combining  numerical  values  where  possible,  we  have 


as(»a) 

at 


21.6 Tll2AC2k3kJ1'2  + 19.44C2*s*,-41/,  - 
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» / 


This  expression  has  been  checked  against  a direct  numerical  integration  of  Eq.  (31)  for 
C - 0.01  and  the  exact  value  of  A . The  results  are  given  in  Table  2.  For  zero  wind  speed, 
Eq.  (33)  becomes 


- 36 «C*r4g-7'V'a  - 207CTf"1/2*|/a  . 


(34) 
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Table  2 — Comparison  of  Energy  Fluxes  Predicted  by  Eqs.  (31)  and  (33) 


*3. 

(cm-1) 

«* 

(cm/s) 

dS(k3)/dt  (Eq.  31) 

(cm4) 

aS(*3)/dt  (Eq.  33) 

(cm4) 

0.11 

16 

2.96  X 10'9  -(1.93  X 10~2)S(k3) 

3.81  X 10'9  -(2.35  X 10'2)S(*3) 

0.11 

60 

1.84  X 10-8  -(3.28  X 10'2)S(k3) 

2.71  X 10‘8  -(4.11  X 10'2)S(k3) 

1.16 

16 

6.87  X 10~8  — (7.06)S(k3) 

9.53  X 10'6  -(7.76)S(*g) 

1.16 

60 

2.73  X 10‘4  -(11.2)S(fe3) 

4.51  X 10~4  - (13.9)S(feg) 

' 

5 


SUMMARY  AND  DISCUSSION 

We  have  assumed  that  < 3.64  cm-1  and  derived  a simplified  form  for  the  gravity- 
capillary wave  interaction.  For  the  zero  wind-speed  case,  this  assumption  is  sufficient  to 
reduce  the  form  of  the  interaction  to  the  simple  analytical  expressions  given  in  Eqs.  (29)  and 
(34).  The  first  term  in  each  expression  is  several  orders  of  magnitude  smaller  than  the 
second  so  that  energy  is  always  drained  from  the  short-gravity  wave. 

To  include  wind-speed  effects  in  the  problem,  the  dispersion  relation  was  modified  as 
given  in  Eq.  (5).  This  procedure  is  by  no  means  rigorously  justified  but  is  an  attempt  to 
include  the  predominant  effect  of  air  flow  in  the  problem,  considering  that  a rigorous  analy- 
sis of  wave-wave  interactions  in  a shear  flow  is  extremely  complex  and  currently  evades 
elucidation. 

Three  additional  approximations  were  necessary  to  simplify  solutions  for  a wind- 
dependent  dispersion  relation.  Two  of  these  are  given  in  Eqs.  (19)  and  (27);  the  third  con- 
sisted of  setting  cos  af  equal  to  one  in  all  dispersion  relations.  Such  dispersion  relations  are 
good  approximations  for  small  angles  where  wave  heights  and  energy  transfers  are  large.  The 
results  of  these  manipulations  are  given  in  Eqs.  (28)  and  (33).  Table  2 shows  that  the  first 
two  terms  in  these  equations  are  very  small.  Furthermore,  the  last  terms  in  the  equations 
are  less  than  about  10%  of  the  third  terms.  Thus,  a rather  good  approximation  to  the  gravity- 

wind  speed  may  be  written 

- -ST.BwgfcjifefS^g.OJS^.O),  (36) 


-92.2wg  *8*3$(*3  )$(*,),  (36) 

where 

3V^/2+  Ut  - s/g/k3  + Ug.  (37) 

To  this  extent,  then,  the  gravity-capillary  wave  interaction  for  short-gravity  waves  may  be 
considered  to  be  an  interaction  between  a short-gravity  wave  and  a capillary  wave  whose 
group  speed  equals  the  gravity -wave  phase  speed. 


capillary  interaction  for  an  arbitrary 

as(»8,0) 

dt 

*S<»  a) 

dt  " 
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